Superconducting proximity effect in interacting quantum dots revealed by shot noise 
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Abstract 

We study the full counting statistics of charge transport through a quantum dot tunnel-coupled to one normal and one 
superconducting lead with a large superconducting gap. As function of the level detuning, there is a crossover from a 
regime with strong superconducting correlations in the quantum dot to a regime in which the proximity effect on the 
quantum dot is suppressed. We analyze the current fluctuations of this crossover in the shot-noise regime. In particular, 
we predict that the full counting statistics changes from Poissonian with charge 2e, typical for Cooper pairs, to Poissonian 
with charge e, when the superconducting proximity effect is present. Thus, the onset of the superconducting proximity 
effect is revealed by the reduction of the Fano factor from 2 to 1. 

Keywords: A. Superconductors, A. Metals, C. Electronic transport, C. Noise, C. Tunneling 
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1. Introduction 

Current fluctuations in a nanoscale conductor yield 
more information than the average current on the pro- 
cesses involved in transport and on the nature of the current- 
carrying excitations. For this reason the study of shot 
noise has attracted a lot of attention both theoretically 
|and experimentally 0, The full current distribution 
function is given in terms of the Full Counting Statistics 
(FCS), P(N,t), defined as the probability that N charges 
have passed through the system in a given time t d 13] • 
The FCS can be conveniently expressed in terms of the 
Cumulant Generating Function (CGF), 



in interacting nano-scale systems @, d 0, @, H E3, EL El 
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S(X) = l n 



^ exp(iN X )P(N,t) 



,N=-oo 



(1) 



from which the current cumulants are obtained by per- 
forming derivatives with respect to the counting field x, 
i.e. «W = (e n /t) d? x S(x)\ x=0 , with e < being the elec- 
tron charge. Several theoretical approaches have been de- 
veloped to determine the CGF both in non-interacting and 



The measurement of high-order moments 
of the current is a highly-demanding experimental task 
0, E, 0, S HIH ■ A breakthrough has been recently 
achieved by measuring in real time single-electron tunnel- 
ing through a quantum dot by means of a nearby quantum 
point contact $H 0, [H, ES ESj ■ The FCS is an important 
tool to st udy transport in hybrid normal-superconducting 
systems [H, 2^, 30, 31, 32, [33| , since in this case it be- 
comes of fundamental importance to identify whether the 
correlations induced by the pairing in the superconductor 
persist in the normal regions. 

Hybrid superconducting-normal structures with quan- 
tum dots offer the possibility to study the interplay of su- 
perconducting correlations, strong Coulomb interaction, 
quantum fluctuations and non-equilibrium physics as a 
function of tunable system parameters. Subgap trans- 
port in these structures is mediated by Andreev reflection 
glliand it has been studied extensively both theoreti 
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[62J,|63|,|64J,|65J,|66|,|67|. 

In this letter, we study the FCS of electrons trans- 
ferred through a single-level quantum dot tunnel coupled 
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to one normal and one superconducting lead with a large 
superconducting gap. The generalized Fano factor Fj = 
K I (~ e) 3-1 /^ 1 ) for Andreev reflection at a normal-super- 
conductor interface is Fj = 2 2 , indicating the transfer of 
Andreev pairs jglj • As we will discuss below, the Fano fac- 
tors remain Fj = 2 J ~ 1 for many transport regimes when a 
quantum dot is placed in between. For strong tunnel cou- 
pling between dot and superconductor, however, the Fano 
factor can be reduced to Fj = 1. In that case, indeed, 
the proximity effect leads to fast coherent oscillations of 
Cooper pairs between dot and superconductor, which are 
stochastically interrupted by injection of single electrons 
from the normal lead. Once the proximity effect is sup- 
pressed either by a large energy detuning of the empty 
and double-occupied dot, or by a large tunnel coupling to 
the normal lead, the Fano factors goes back to Fj = 2 3 ~ l . 
This counter-intuitive result shows that the superconduct- 
ing proximity effect on the dot can be revealed by a mea- 
surement of current fluctuations, rather than the current. 



2. Theory 



the limit A — > oo, in which the superconducting lead's de- 
grees of freedom can be integrated out. This results in the 
effective Hamiltonian that for the dot plus the supercon- 



ducting lead [41 



H, 



dot 



t (44 



(5) 



The tunnel coupling to the superconductor mixes the states 
|0) and \D) to form new eigenstates |±) = ^jaA T 5^|0)-F 

-±=Jl±^-\D) with energies 6/2 ± e A . Here 6 = 2e + U 



s/2 V * "~ 2e 

is the detuning between empty and doubly-occupied state 
and 2e A = \/6 2 + T| the splitting between the |+) and 
the |— ) states. The excitation energies of the dot define 



the Andreev addition energies 



E An , n = -f , U/2 + 1 e A , 



(0) 



with 7', 7 = ±1, whose position can be identified by mea- 
suring the charge or the current as a function of the elec- 
trochemical potential [i [Hl|, as long as Yjq <C ksT <^ e A . 



We consider the Hamiltonian H = H do t+J2 n =N,s( H v+ 3. Calculation 



Htunn,ri), where r\ = TV, S 1 labels the leads. The dot is 
described by the Anderson impurity model with a spin- 
degenerate level £ and on-site Coulomb repulsion U 



(2) 



where d a {d] J ) is the annihilation(crcation) operator for an 
electron with spin a =t, 1 and n a = d\d a . The dot eigen- 
states are denoted as |0) (empty), \a) = d£|0) (singly oc- 
cupied with spin er), and \D) = d|d||0) (doubly occupied). 
The Hamiltonian of lead 7/ = TV, S reads 

H n e^c^Crjkc ~ S, hS ^2 (c n -kic v kt + H.c.) , (3) 

ka k 

with the lead-electron operators c n ka and c\ ka , and 5„ s = 



"7]ka ' 

1 for 77 = S and for 77 = N . Superconductivity in lead S 
is described by means of the mean-field BCS Hamiltonian 
with a pair potential A = fl l Sfc( c ifct c i-fei)' wncrG 9 de- 
notes the interaction strength. Without loss of generality, 
A has been taken real and positive. We choose the elec- 
trochemical potential of the superconductor as a reference 
/js =0. The normal lead has electrochemical potential 
fiN = f/,. The tunnel coupling between dot and leads is 
modeled by the tunneling Hamiltonians: 



H 



tunn,?7 



V, 



ka 



H.c. 



(4) 



We assume that the normal-state density of states v n of 
lead 77 is constant in the energy window relevant for trans- 
port and define the tunnel-coupling strength as r ?) 
Since we are interested in sub-gap transport, we consider 



2-Kv n \V n \ 



We trace out the normal lead's degrees of freedom in 
order to obtain an effective description of the dot in terms 
of its reduced density matrix. For Tn ^ e A , the reduced 
density matrix is diagonal in the cigenbasis of the effec- 
tive Hamiltonian ([5]). The vector of the diagonal matrix 
elements, the occupation probabilities, at time t under the 
condition that TV electrons have passed through the bar- 
rier between the normal lead and the dot, p(TV, i), obeys 
the generalized master equation (l2l | 

jP(N,t)= / dt'W(N-N',t-t')p(N',t'), (7) 

where the Kernel W(TV— TV', t—t') describes transitions in- 
volving the transfer of TV— TV' electrons. The counting field 
is introduced through the Fourier transform W x (< — t') = 
J2n W(TV, t — t') exp(ixTV). In the paper we consider only 
the weak-coupling regime Tn -C ksT and we perform a 
perturbation expansion to first order in Tat. To this or- 
der of perturbation theory non-Markovian corrections play 
no role 12j and we need only the zero-frequency Laplace 
transform of the Kernel: W x = J Q ~W x (r)dT. The long- 
time behavior is determined by the eigenvalue of W x with 
the smallest absolute value of the real part, A(x) [HI, HI]. 
In this regime the CGF in first order in r^v is given by 
S = tX*- 1 ' (x), where t is the measurement time and the su- 
perscript explicitly indicates first order in Yjy. A general 
diagrammatic theory for computing the Kernel W x exists 
(5l| . but in the case under consideration the rates W^, 
for transitions from the dot state 7// to the state tp can be 
' computed straightforwardly by means of Fermi golden rule 
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Figure 1: Schematic description of the rates contributing to transport 
in the high-bias regime, a) On resonance, all rates are the same and 
all give rise to the transfer of one electron from the normal lead to 
the dot. The states |±) are coherent superposition of |0) and \D) 
with equal weight: |±) = (|0) =F -D))- b) Off resonance there are 
fast processes (solid lines) and slow processes (dashed lines). For 
<5 > r s , the dot will be mostly stuck in the state |+> ~ ^f|0) + \D). 
Transport occurs along the red cycle (lighter gray) and it consists of 
uncorrelated tunneling of pairs of electrons. 



introducing a factor exp(— i\) [cxp(ix)] for each electron 
leaving [entering] the normal lead. The rates read: 



7=±1 v 



- 7lT7 )e-^ (8) 



7=±1 



A, 7, ±7 ) 



-»7X 



(9) 



where f + (j-o) is the Fermi function of the normal lead 
with electrochemical potential [i, and /"(w) = 1 — f + (uj). 
The diagonal elements of the Kernel do not contain the 
counting field and are given by W+,+ = — 2W* + | x= o, 
= -2W*_| X=0) and W a ,„ = -(W* j(y + W*J)\ X=Q . 



4. Results 



r| for which the 



We discuss the regime \S\ < \/TP 
Andrecv bound states obey the inequality Ea,+,+ > Ea,+.- 
Ea,-,+ > Ea.-.-- In the following we consider only the 
low-temperature limit fc^T <C ca and the case of positive 
/i. The results for negative fi can be obtained by per- 
forming the substitutions 6 — 8 and x ~ X- In the 
high-bias regime, fi > E A ,+ t +, the CGF reads 



5; 



hb 



tr 



N 



r|e- 2 *x 



8 2 



(10) 



On resonance, |(5| <C Ts, this expression simplifies to the 
CGF corresponding to Poissonian tunneling with strength 
Fjv of electrons through a barrier 



S- 



hb.on 



tT N (exp(-ix) - 1) , 



(11) 



Andreev current hK^/(eT N ) 
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Figure 2: Current in the normal lead in units of eTpf/h as a function 
of level position e and of chemical potential fi for = 0.2 U, tjT = 
0.01 U and Tat Tg. The dashed lines map the Andreev bounds 
states as fi = E a 



while off resonance, |<5| ^> Ts, we obtain 

r 2 

S hh , oS =tr N ^(exp(-2ix)-l), (12) 

which corresponds to a Poissonian transfer of Andreev 
pairs. In between the two limits, the FCS is non-Poissonian. 
We now proceed to clarify the physical picture that leads 
to this somewhat counter-intuitive result. 
On resonance (\S\ <C Tjv): In the large- bias regime, the 
tunneling rates are all equal and they are given by 



(T N /2) expHx). 



(13) 



The situation is pictorially shown in panel a) of figure [TJ 
The stationary probability of being in the |+), |— ) and 
\a) states are the same and equal to 1/4. Suppose that 
> at a given time the dot is in the singly-occupied state \<j). 
On average, after an interval of time of 2h/T n an electron 
from the normal lead enters the dot. The resulting doubly- 
occupied state, \D) = (|— ) — |+))/V2, undergoes fast co- 
herent oscillations to |0) and back, describing tunneling of 
an electron pair forth and back between quantum dot and 
superconductor. This happens well before than another 
electron enters in the dot after a time 2h/T^. Seen from 
the normal conductor, single electrons hop on the dot in- 
dependently with a rate Fat/2. 
Off resonance (\8\ 3> Tat): For simplicity we discuss only 
the 8 > case. The rates are given by 



w* >t 



r Ar r|/(4<5 2 )exp(- ?x ) 

T N exp(-ix) , 



and the dot will be with a high probability in the state |+) . 
Panel b) of figured] elucidates the situation. Notice that in 
this regime the state |+) corresponds mostly to \D) with 
only a small admixture of |0). The dot first needs to inject 
a Cooper pair in the superconductor via a cotunneling pro- 
cess (factor r|/(5 2 ) and only then an electron can tunnel 
in the dot from the normal side (factor T^.) This com- 
bined process is the bottleneck for transport and leads to a 
suppression of the rate W* + , which goes as TnT^/S 2 . On 
the contrary, when the dot is in \a), the rate for tunneling 
back to |+) is much higher oc Tn and the dot will return 
quickly to the state |+). Seen from the normal side, pairs 
of electrons arc injected with an average time separation 
of 4M 2 /(rjvr|). The average separation of the electrons 
within a pair is much shorter, namely h/T^. 

In the intermediate-bias regime, £U.+,- < M < Ea,+,+, 
the situation is more complex. The CGF reads 



Sib — t 



i_N_ 



S-6e A + J [35 + 2e A ) 2 + 8T 2 s e 



-2i X 



(14) 



The Poissonian one-electron transport regime is now re- 
alized for 5 = — ^T,s where the CGF becomes = 
t^f- (exp(— ix) — !)• Off resonance, we recover the Pois- 
sonian tunneling of Andreev pairs 



Si 



ib,off 



tr 



i 



N 



S 2 3 + sgn 5 



(exp(-2i X ) - 1) 



(15) 



This behavior can be understood looking at the rates in 
a similar fashion as for the high-bias case. Finally, in the 
low-bias regime, < /i < Ea,+.-, transport is suppressed 
by Coulomb blockade. In figure [5] we show the current 
that, as already discussed in reference [Hl|], exhibit a de- 
pendence on the Andreev addition energies (dashed lines) 
as a function of e and /i. 

In figure [3l for fixed bias, we show the transition from 
Fano factors Fj sa 2^ 1 to Fj « 1 for F 2 (black) and F 3 (red) 
as a function of e. Off resonance, we have Fj « 2- 7 " 1 , in- 
dicating that electrons are transfered as Andreev pairs. 
On resonance, the superconducting proximity effect leads 
to fast coherent oscillations of a Cooper pair between dot 
and superconductor. Electrons are injected from the nor- 
mal lead via uncorrelated single-electron tunneling events, 




Figure 3: Fano factor F2 (black lines) and normalized skewness F3 
(red lines) as a function of level position for fi = 1.5 f/ (solid line) 
and fA = 0.5 U (dot-dashed line). The other parameters are the same 
as in figure [2] 



extreme limit of Tn 3> £a , the superconducting proxim- 
ity effect is fully suppressed, and the picture of uncorre- 
lated electron-pair tunneling remains valid for all values of 
6/Tg. We explicitly checked the cases r^r 3> \S\ 3> and 
Tjv 3> Tg ^> \6\ in the high-bias regime (by means of the 
diagrammatic theory (Hlj and including matrix elements 
of W that are off-diagonal in the |±) basis) to find the 
CGF 



S = t- 



2T 



(exp(-2 lX ) - 1) , 



(16) 



that corresponds again to Poissonian transfer of Andreev 
pairs with the rate r|/2Fjv- For intermediate values of 
Tjv, representing finite but not complete dephasing, we 
expect an intermediate value of the Fano factor between 1 
and 2 on resonance. In that sence, the measurement of the 
Fano factor may serve as a tool to determine the presence 
and strength of dephasing. 

5. Conclusions 

We have computed the CGF for transport through a 
quantum dot coupled to one normal and one superconduct- 
ing lead in the large-gap limit. In the high-bias and inter- 
mediate bias regimes, there is a crossover from generalized 
Fano factors Fj = 2 5 to Fj = 1, which clearly identifies 



Fi ~ 1. In the intermediate-bias regime (dot-dashed lines , u j, ,■ ■ ■, n(t , • ,1 , , , 

, . . . the superconducting proximity effect m tne quantum dot. 



in figure [3]) we see a similar behavior but with a shifted 
position of the minimum. 

In the analysis presented above no dephasing giving 
rise to a suppression of the superconducting proximity ef- 
fect is included. In fact, in the limit of weak tunnel cou- 
pling to the normal lead, <C £a, as assumed above, the 
Andreev bound states are sharp resonances without broad- 
ening. The regime of a quantum dot more strongly coupled 
to a superconductor than to a normal lead, i.e. T$ ^> T^r, 
has already been realised experimentally by Deacon et al. 
[H, |67|]. The effect of dephasing can be studied by in- 
creasing the tunnel coupling to the normal lead. In the 
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